Spin gradient demagnetization cooling of ultracold atoms 
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We demonstrate a new cooling method in which a time- varying magnetic field gradient is applied 
to an ultracold spin mixture. This enables preparation of isolated spin distributions at positive and 
negative effective spin temperatures of ±50 picokelvin. The spin system can also be used to cool 
other degrees of freedom, and we have used this coupling to cool an apparently equilibrated Mott 
insulator of rubidium atoms to 350 picokelvin. These are the lowest temperatures ever measured in 
any system. The entropy of the spin mixture is in the regime where magnetic ordering is expected. 
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Attainment of lower temperatures has often enabled 
discovery of new phenomena, from superconductivity to 
Bose-Einstein condensation. Currently, there is much in- 
terest in the possibility of observing correlated magnetic 
quantum phases in lattice-trapped ultracold atoms [TH3] • 
The relevant critical temperatures are on the order of 200 
picokelvin, lower than any previously achieved tempera- 
ture jU |5j . Realization of this temperature scale requires 
the development of new methods of refrigeration which 
can be applied to ultracold atoms. Many such techniques 
have been proposed [6T-I13] but await experimental real- 
ization. The cooling method we demonstrate here opens 
up a previously inaccessible temperature regime and pro- 
vides a realistic path to the observation of magnetic quan- 
tum phase transitions in optical lattices. 

The new cooling method is applied to an optically 
trapped cloud of cold atoms in a mixture of two internal 
states with different magnetic moments |14j . Application 
of a strong magnetic field gradient results in almost com- 
plete spatial segregation of the two spin components. The 
"mixed region" of spins between the pure-spin domains 
has a width which is proportional to the temperature T 
and inversely proportional to the applied gradient. Re- 
ducing the gradient mixes the two components, and due 
to the mixing entropy the temperature is dramatically 
reduced. Since Mott insulators can be prepared with an 
entropy per particle much lower than fcs In 2, our scheme 
is a practical way of creating a low entropy mixture in the 
regime where magnetic ordering is expected. This cooling 
scheme introduces several new concepts. It implements 
demagnetization cooling with spin transport instead of 
spin flips, allows decoupling of the spin and kinetic tem- 
peratures, and enables the realization of negative spin 
temperatures. Long tunneling times in the lattice allow 
two different implementations of our cooling scheme. 

When the gradient is changed faster than the spin re- 
laxation rate, the spin system is effectively isolated from 
all other degrees of freedom, and very low spin tempera- 
tures can be achieved. Contrastingly, when the gradient 
is changed slowly enough, the spin system is fully equi- 
librated and can absorb entropy from other degrees of 
freedom, cooling the whole sample. Reduction of the 
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FIG. 1: Two different cooling protocols realizing the cases 
of isolated and equilibrated spin systems, a: Experimen- 
tal "phase diagram" of lattice depth vs. applied gradient. 
Dashed lines show two different paths which connect the high- 
gradient superfluid state and the low-gradient Mott insulat- 
ing state, b: and c: show the lattice depth (solid line) and 
gradient strength (dashed line) versus time for the two cases 
(equilibrated spins and isolated spins, respectively) in panel 
a. The shape of the lattice rampup is designed to ensure 
maximum equilibration. 



gradient after (before) the optical lattice has been raised 
realizes the regime of isolated (equilibrated) spins. 

First, we discuss isolated spins, of which atoms in a 
Mott insulating state are an almost ideal realization. 
Spin distributions relax by two atoms exchanging sites 
through a second order tunneling process. The time scale 
of this relaxation is typically one second, and the gradi- 
ent can easily be varied much faster. The equilibrium 
spin distribution depends only on the ratio of the ap- 
plied gradient V|B| and temperature T. When V|B| is 
changed the effective temperature of the decoupled spin 
degrees of freedom (spin temperature) is rescaled pro- 
portionally. This enables the realization of spin distribu- 
tions with a very low positive (or, if the sign of the gra- 
dient is changed, negative) effective temperature. Nega- 
tive temperatures can only occur for systems with an up- 
per bound on the energy, and correspond to an inverted 
Boltzmann distribution with the largest population in 
the highest-energy state [15] , 

Figure [2] shows the results of experiments on isolated 
spins. Fits to data on equilibrated spins indicate an ini- 
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tial temperature of 6.3 nanokelvin (see Figs. [3] and [4j. 
While reduction of the gradient by a factor of 1000 would 
be expected to reduce the effective temperature of the 
isolated spins to 6.3 picokelvin, our finite optical resolu- 
tion only allows us to assert an upper bound of 50±20 
picokelvin. We have held the spins for up to several 
seconds in the lattice. No heating is observed during 
a one-second hold, but after three seconds the 50 pi- 
cokelvin distribution is observed to heat to about 70 
picoKelvin. Similarly, at small negative gradients, we 
observe a negative-temperature distribution with a tem- 
perature closer to zero than -50±20 picokelvin. Since the 
total energy is monotonic in — 1/T, these are among the 
most extreme thermodynamic states ever measured in 
any system [TH1 E] . Very low condensate release energies 
(not temperatures) have been observed previously [15] , 
Note that the spin temperatures we report are much 
lower than those attainable by magnetic trapping or opti- 
cally pumping a system into a single spin state: even for a 
fractional population of 10~ 6 in other spin states caused 
by imperfect pumping or spin-flip collisions |19| , the spin 
temperature in a bias field of 100 mG is 500 nanokelvin, 
assuming a magnetic moment of one Bohr magneton. In 
our experiment, the energy scale is set by the product of 
the Bohr magneton with the magnetic field gradient and 
the lattice spacing, and it is the relative ease of achieving 
very small values of this parameter (corresponding to \xb 
times a field of less than one microgauss) which allows us 
to reach such low spin temperatures. 
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FIG. 2: Preparation of low positive and negative spin tem- 
peratures. Measured spin temperature versus final gradient, 
for the case of isolated spins. Error bars are statistical. 

The mixed-spin region comprises a spectrum of spin 
excitations, the energy of which can be tuned by adjust- 
ing the strength of the gradient J20|. A spin excitation 
consists of a spin-up atom swapping places with a spin- 
down atom on the other side of the zero-temperature spin 
domain boundary. Reduction of the gradient reduces the 
energy of spin excitations. In the regime of equilibrated 
spins, this causes entropy to flow into the mixed-spin re- 
gion from other degrees of freedom. This lowers the tem- 
perature of the whole system in a manner locally anal- 
ogous to standard single-shot adiabatic demagnetization 



refrigeration [5TJ [52]. The kinetic excitations of a trapped 
Mott insulator are particle- hole excitations [4[ [23] . On 
a microscopic level, particle-hole excitations can couple 
directly to spin excitations if a particle and a hole on 
opposite sides of the spin domain boundary annihilate. 
Energy may be absorbed by quasiparticles in the first 
band or by the spin excitations themselves [21] • In the 
superfluid phase and during lattice rampup, the kinetic 
excitations are different and entropy transfer is expected 
to occur faster and in more complex ways. Experimen- 
tally, it is easier to maintain adiabaticity if most of the 
path along which the gradient is changed is in a regime 
of fast relaxation times (e.g. the lower path in Fig. la). 
Thermodynamically, this adiabatic cooling method is a 
redistribution of entropy from the kinetic degrees of free- 
dom to the entropy which results from partially mixing 
the two initially segregated spin domains. 

It is easily possible for the mixed region to absorb 
nearly all of the entropy of the system. In a one- 
component harmonically trapped Mott insulator which 
is at a temperature low enough for the particle-hole ap- 
proximation to hold, the approximate total entropy is 
fee ln(2) times the volume of the "shells" between the 
Mott plateaux [3]. The maximum entropy of the the 
mixed region is realized when, at low gradient, it is broad- 
ened to a substantial fraction of the total size. In that sit- 
uation the entropy per site approaches k-Q ln(n+l), where 
n is the local number of indistinguishable bosons per site 
(in our samples, n varies across the trap between 1 and 3). 
The maximum entropy of the spin system is thus larger 
than the entropy of the kinetic degrees of freedom. Thus, 
substantial cooling of the system can be achieved with a 
single gradient demagnetization ramp. We have made a 
more quantitative analysis of spin gradient demagneti- 
zation by calculating entropy-versus-temperature curves 
of our system in various field gradients [25]. The results 
confirm the qualitative argument above and show that 
spin gradient demagnetization cooling is capable of cool- 
ing well below the expected Curie temperature. 

Figures [3] and [4] show the results of spin gradient de- 
magnetization cooling experiments. As the gradient is re- 
duced, the width of the domain wall increases (see Fig.|3|, 
indicating transfer of entropy from the kinetic degrees 
of freedom to the spins. The width increases much less 
steeply than would be expected for an isothermal sample, 
implying cooling. The observed domain wall width can 
be converted to a temperature using spin gradient ther- 
mometry. The measured temperature falls rapidly as the 
gradient is lowered (see Fig. [1]). The lowest measured 
temperature is 350±50 picokelvin, making this the cold- 
est implementation of the atomic Mott insulator. An im- 
portant goal for future experiments will be to develop an 
additional method of thermometry to measure the tem- 
perature of the kinetic degrees of freedom at this very 
low temperature scale (for progress towards this goal, see 
Refs p?6Ti28j ). This temperature is colder than the low- 
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FIG. 3: Entropy transfer from other degrees of freedom to the 
spin system. Circles represent measured width of the mixed 
region vs. final magnetic field gradient, for the "equilibrated 
spins" protocol (see Fig. [I]). Error bars are statistical. The 
dashed line represents the expected behaviour assuming no 
cooling. The dash-dotted line shows the minimum resolvable 
width. The solid curve is the theoretical prediction, assuming 
an initial temperature of 6.3 nanokelvin and including the 
effects of optical resolution (see Ref. [25] for details). Insets 
show spin images at the indicated points, the corresponding 
vertically integrated spin profiles (squares), and the fit to the 
expected form of a tanh function times the overall density 
distribution (solid line; see Ref. (20] f° r details). Axis units 
and grayscales in the two insets are abitrary but identical. 



est temperature ever measured in an equilibrated kinetic 
system [2!5] , and it is within a factor of 2 of the expected 
magnetic ordering temperature [5J. 

Theoretical curves in Figs. |3] and ( see also refer- 
ence [25]) show reasonable agreement with the data. 
These curves were fitted to the measured temperatures 
using only one free parameter: the initial temperature 
at the maximum gradient. The initial temperature in- 
ferred from this fit is 6.3 nanokelvin. In our earlier work 
on thermometry [20] . where the lowest measured tem- 
perature was 1 nanokelvin, some adiabatic demagnetiza- 
tion cooling may have occurred during the preparation of 
the system. The flattening-out observed in the measured 
data at low gradients could be a signal that all accessible 
entropy has been pumped into the spin system. 

There are both practical and theoretical limits on the 
temperatures which can be attained with spin gradient 
demagnetization cooling. In traditional magnetic refrig- 
eration experiments, the minimum temperature is often 
set by the minimum achievable magnetic field or the pres- 
ence of internal fields in the refrigerant [30]. Analogues 
of both these limits are relevant here. Practically, the ra- 
tio between the highest and lowest attainable magnetic 
field gradients is an upper bound on the ratio between 
the initial and final temperatures. In our experiment, the 
maximum value of V|Bj|/V|B/| is about 1000 (limited 
mainly by the accuracy of determining the zero crossing 
of the gradient), which would give a minimum tempera- 
ture below 10 picokelvin. Another limit stems from the 
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FIG. 4: Spin gradient demagnetization cooling. Circles repre- 
sent measured temperature vs. final magnetic field gradient, 
for the "equilibrated spins" protocol (see Fig. [T]) . Error bars 
are statistical. These measurements are the same as those 
shown in Fig. [3] The dashed line follows the isothermal tra- 
jectory and the dash-dotted line shows the resolution limit. 
The solid line is the theoretical prediction, assuming an initial 
temperature of 6.3 nanokelvin and including the effects of op- 
tical resolution. The dotted line is the theoretical prediction 
without the effects of optical resolution. 



small difference between the inter-spin interaction energy 
U-fi and the mean of the intra-spin interaction energies 
(E/ T f + Un)/2. For 87 Rb, this limit is not expected to 
preclude cooling below the expected magnetic ordering 
temperature [25] . The fundamental limit is set by the 
total entropy of the system. A single Mott insulator has 
an entropy per particle much less than fee In 2. At high 
gradient, the spin entropy is negligible. Lowering the gra- 
dient provides a controlled way of creating a completely 
mixed two component system at the same entropy which 
is low enough for spin ordering to occur. This technique 
thus provides a specific and realistic method of realizing 
magnetic phase transitions in optical lattices. 

Our scheme differs from single-shot adiabatic demag- 
netization refrigeration (including that demonstrated in 
a gas of chromium atoms [3T]) in that the magnetic field 
is replaced by a magnetic field gradient, and spin flip col- 
lisions by spin transport. The chromium scheme cannot 
be applied to alkali atoms due to the much slower spin-flip 
rates, and extending it from microkelvins to picokelvins 
would require sub-microgauss magnetic field control. In 
previous work [20], we suggested that adiabatic reduc- 
tion of the gradient could be used for cooling, and some 
aspects of this proposal have recently been theoretically 
addressed and verified [32 . 

The concept behind spin gradient demagnetization 
cooling is compelling; if adiabaticity can be maintained, 
then strong cooling in a lattice will occur. Our experi- 
mental implementation was designed to allow the system 
to equilibrate as much as possible at low lattice depths. 
We have tested reversibility by replacing the single gra- 
dient ramp by a sequence of ramp down, ramp up, ramp 
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down. This led to no detectable difference in the final 
measured temperature. This indicates that the gradient 
ramps are adiabatic. Equilibration in the Mott insula- 
tor is more difficult to demonstrate, although the pre- 
viously demonstrated agreement between spin gradient 
thermometry and cloud size thermometry at high tem- 
peratures [20 indicates that the kinetic and spin degrees 
of freedom are equilibrated in that regime. The fact that 
the spin distribution fits well to the form expected of an 
equilibrated spin system is also evidence for equilibration, 
as is the one-parameter fit to our theoretical predictions 
(which assume adiabaticity) . However, if the lattice is 
deepened, then lowered to zero, then raised again, heat- 
ing is generally observed. Thus, we cannot rule out the 
existence of long-lived metastable excitations in the Mott 
insulating state which do not couple to the spin degrees 
of freedom and thus do not influence our temperature 
measurement. Other experiments have seen evidence of 
long equilibration times in the Mott insulator [211 l33j . 
For quicker equilibration, spin gradient demagnetization 
cooling could be implemented with lighter atomic species 
(e.g. 7 Li or 4 He*) and/or shorter period optical lattices. 

The cooling technique using isolated spins presented 
here has achieved spin temperatures and entropies well 
below the critical values for magnetic ordering, and spin 
gradient demagnetization cooling of equilibrated spins 
has cooled to a point within reach of the critical values. 
This work thus opens a realistic path towards observa- 
tion of superexchange-driven phase transitions in opti- 
cal lattices, and extends the potential of ultracold atoms 
trapped in optical lattices to be used as flexible quantum 
simulators of strongly interacting many-body systems. 
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